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Abstract 

The density matrix is obtained for short-wavelength photons passing through a 
medium. For sufficiently large thicknesses of a target this matrix is expressed by the 
forward scattering amplitude of a photon. For multi-GeV photons scattering via 
virtual e + e~ pairs is the most important process. Its amplitude in a crystal is cal- 
culated. The results obtained are used to find the optimal thicknesses, orientations 
and types of crystals for the circular - to - linear polarization conversion process. 



Introduction 



High-energy circularly polarized photons, if available, would help in solving several im- 
portant problems in high energy physics like that of spin-crisis (see e.g.Q and literature 
cited there). If we start with linearly polarized multi-GeV photons, then they can be 
converted into circularly polarized ones provided that the suggestion of Cabibbo and col- 
laborators to use a specially chosen crystalline plate for this purpose is true. This idea 
is now under experimental investigation at CERN within NA59 project, where the linearly 
polarized photons are also produced by means of a crystalline target from unpolarized 
electrons passing through it. From theoretical point of view, the conversion process itself 
is very interesting since for multi-GeV photons it is caused by the polarization of vacuum 
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in the presence of the periodical electric field of a crystal. This phenomenon has not been 
investigated so far. 

To develop a description of the photon propagation, we start in Sec. 2 from the Maxwell 
equations taking into account a current induced in a medium by the incident wave in a 
rather general form. When the wavelength of a photon is much shorter than any other 
characteristic distance scale of the problem, the parameters of a wave packet change 
slowly while it propagates in a medium. Using this fact, we obtain the solution to the 
Maxwell equations in the short-wavelength approximation, which is the analog of the 
eikonal approximation in the fast particle scattering theory. The approximation used 
becomes valid starting with a relatively small photon energy u > m , where m is the 
electron mass ( we use h = c = 1 units). Comparing our solution with that of || 
obtained for an amorphous medium, we find that a matrix, describing a change of the 
polarization and intensity of a wave packet is deeply connected with the forward scattering 
amplitude. At least, it is the case for large thicknesses, when this change becomes really 
prominent. In the multi-GeV energy region of interest, the photon scattering via virtual 
e + e~ pairs is the only process relevant to the problem. Its amplitude for a separate atom 
is well known ( see [|J for the forward scattering amplitude ). In crystals, along with this 
incoherent ( amorphous like ) contribution to the amplitude, the coherent one, caused by 
the periodicity of a lattice is present. The latter is calculated in Sec. 3 by means of the 
so called quasiclassical operator method. The details of this method along with many 
applications can be found in || . We use the results obtained to consider the circular - to - 
linear polarization conversion process. Optimal ( according to the criterion formulated in 
H ) thicknesses and orientations are found for diamond, silicon and germanium crystals 
for photon energy w ~ 100 GeV. 



2 Propagation of short-wavelength photons 



For the electric field E of a wave we obtain from the Maxwell equations 



d 2 
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where the operator 



( - V 2 )E t ( X ) = 47T J (Px'RfrnWXmjfax^Ejtf), (1) 
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When deriving eq.(^), the following relationship of the electric E and the induced D fields 
was used: 

Di{x) = Ei{x) + 4vr J d 4 x' X ij(x, x) E 3 (x') , (2) 

as well as the condition dD^x) / dxi = 0. By definition, differentiating the integral in 
eq.(Q) over t , we obtain j(x) being a density of the induced current. When the function 
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Xij(x,x') depends only on the difference t — t' , the equation for the vector potential A 
is essentially the same as eq.([l|) for the electric field. 

In the vacuum the solution to eq.(0), satisfying the condition dEi(x) / dxi = , reads 

E (x) = / dkg(k - k )e ± (k)e- ikx , (3) 



where x and k are 4-vectors: x = (t, r), k = {u , k) with k = uv ,u 2 = 1 , and e ± (k) 
being an arbitrary vector, perpendicular to k. If we introduce 5^ = Sij — Vi Uj , then 
e^(k) = S^ej. We assume that the function g in eq.(|D vanishes except of the narrow 
region where | k — k | <C | ko, | . So, the incident wave is a wave packet propagating 
along ko . Let it encounter a medium at some boundary . We try to satisfy eq.(|l]) in the 
medium using the field E in the form of 

Ei{x) = J dkg(k - ko)F y (r, k) ef(k) e~ ikx , (4) 

with Fij(r , k) = dij on the boundary. Assuming that uj is sufficiently large, we expect 
that the function i^-(r , k) varies very slowly with respect to r . Substituting eq.(^) into 
eq.(||) and neglecting the term V 2 F ( keeping only kd/drF), we obtain 

= / dkg{k - k )\2ie- tkx k i mj [ r ' k) + 

(5) 

+ 4ti J d 4 x' Ri m {x) Xm i(x, x') F y (r' , k) e~ ikx '} ef(k). 
The integral / d A x' in eq.(^J) can be easily taken if we use the Fourier-transforms: 

(6) 

^(r',k) = J ^ e ^F y (q,k). 

The most general form of the function Xmi{ki, k%) in crystals reads 

Xmi(kuh) = (2vr) 4 ^2c ml (Q,k l )5(k 2 - k x - Q), (7) 
Q 

where Q = (0, Q) is a reciprocal lattice vector. Using this representation, we obtain from 
eq.(|) 



= J dkg(k - k Q )e- lkx {2ik l dF ^ r r ,k) + 

+ Arc^SL J ^0-^c ml {Q,k + q- Q)e^ F tj (q , k) e^} ef(k). 



(8) 



The second argument of the function c m i in eq.(§) is k\ = k + q — Q = (lu, k + q — Q). 
This function depends actually on uj\ — uj and ni = ki / | ki | « v + s -1 , where 
s = (q — Q)/u. We have already used the fact that | s | <C 1 neglecting it as compared to 
unity. In particularly, that is why we have 8^ in eq.(|]) instead of S im — k Vl k lm luj 2 . We 
shall see below that owing to this replacement the longitudinal components of Fy(r , k) 
do not appear at any depth if they are absent on the boundary. However, when v is 
almost parallel to some crystal axis, the transverse ( with respect to this axis ) component 
of i> should be compared with | s | . We can rule out q — Q from the second argument 
of the function c m i when the angle of incidence $0 with respect to this axis ( transverse 
component of v ) is sufficiently large: $0 S> | s | . In what follows, we assume this 
condition to be fulfilled, bearing in mind that really it does not lead to serious limitations 
as a typical magnitude of | s | at u ~ 1 GeV is 10~ 6 . The integral over q turns 
into Fij(r , k) according to the definition eq.(|6|) , when we change c m ;(Q , k + q — Q) =^ 
c m i(Q,k) in eq.(|8p. This means that the integral over r' in eq.(|JJ) converges in the 
domain | r — r' | <C R where R is a characteristic scale for the noticeable change of the 
function F^(r'). So that, we could substitute Fj(r') =>- Fij(r) yet in eq.(^). Now eq.(|) 
goes over into 

= J dkg(k - k )e- i ^{2ik l dFi f r ^ ) + 

(9) 

+ Ancu^t £ c ml (Q , k)e~ iQr F y (r , k) } ef(k). 
Q 

Here we are interested in the transverse ( with respect to u ) tensor Fij . Let us an- 
alyze, whether such a form of the sought tensor is consistent with eq.@ and bound- 
ary conditions. According to eq.(P), the left longitudinal components of this tensor 
Fj = UiFij(r , k) are independent of the penetration depth. If we suppose that Fj 1 = 
on the boundary, then Fy = 5i n F n j = S^F n j at any depth. In particularly, this im- 
plies that only transverse components of the tensor c m i(Q , k) are present in eq.@. The 
right longitudinal components of the sought tensor do not enter eq. (|]) owing to the factor 
ej"(k) . So, all the tensors below are self-consistently assumed to be transverse. They can 
be presented as two-dimensional matrices. Remembering this, we obtain from eq.(^) 

1 

F(r , k) = exp { 2mujz £ c(Q , k)e' iQ P f dx e~ ixzQ " } , (10) 

Q 

where p = r — zv , z = ur,Q\\ = vQ . For Q\\ = the integral in eq. ([lOD equals 
unity, while for Qu ^ it is of the order of {Q\\z)~ l . Noticeable effects appear when 
the main ( Q\\ = ) term in the phase of eq. (|10[) is of the order of unity. As shown in 
the next Section, it happens for z of several centimeters when the contribution of terms 
with Q|| 7^ to this phase ~ (Q^z)^ 1 1 and can be neglected. Since matrices 
Eq c m ;(Q , /c)e~ lQr are in general non-commutative at different z , the exact solution to 
eq.(|) is the z -ordered exponential function. However, corrections to eq.(|lOD due to z 
-ordering are connected with terms Q\\ 7^ being as small ( (Q^z)^ 1 ) as terms already 
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neglected. It must be remembered also that ( for the coherent yield ) Q is a discrete 
vector and that large values of | Q | do not contribute to the sum in eq. flToD . These 
means that for sufficiently thick crystals we should keep only terms with Q» = in the 
sum of eq. (|10|) : 

F(r , k) = exp { 2niuz £ c(Q± , k)e~ iCl± P} . (11) 

Qx 

The condition Qn = implies that Qx = too, if u is not aligned on some crystal 
plane. More precisely, let ip be the angle of v with respect to the plane. For ip ^> 
4>q = d p i/z , with d p i being the inter-planar distance, we can retain the single term 
c(Q = 0, k) of the sum in eq. (|iT|) . The angle ipo is extremely small, even for very thin 
films. The integration over v in eq.([|) should smear out a dependence of the general 
result on p. We are not interested here whether such a dependence is observable or not. 
In what follows, the sum J2q ± c(Q± , k)e~ tCi± P in eq . (|TT|) will be replaced by the matrix 
c(Q = , k) = c(0,k) which coincides with the average of this sum over p. The matrix 
c(0 , k) has no singularities and is practically constant in the region xjj < ip Q . So that, 
any directions are allowed for u . 

Keeping in the sum J2q ± c (Q_l > k)e~ lCl± P the only term with Qj_ = , we disre- 
gard possible small-angle photon scattering. In this approximation there is no exchange 
between fractions of the wave-packet having different v ,i.e., these photons propagate in- 
dependently. So that, for given direction, we can go over to the conventional idealization 
of a monochromatic plane wave. . Then the electric field is obtained from eq.flTTD and 
eq.([|) with the change in the latter g(k — k ) => 5(k — k ) and then k k 

E(x) = e-^e i ^"e ± (k) , (12) 

where the two-dimensional vectors E(x) and e _L (k) are correspondingly the electric field 
in a medium and the polarization vector of the incident wave. The quantity n = 5^ + 
27rc(0 , k)ij is a complex matrix representing, by definition, the index of refraction. It 
depends on u and ( in crystals ) on v . We introduce a matrix II = 4tiuj 2 c(0 , k) , then 
n — I + U/(2u 2 ) . The photon density matrix at depth L reads 



p(L) = expj^jp(0)exp j-zL^j , (13) 



where p(0) is the density matrix of the incident photon. 

Comparing eq. (|T2"D with the corresponding result of || obtained for an amorphous 
medium by the direct summation of fields coming from individual scatterers, we find that 
in this case II = 4nf(0) N . In this formula f(0) is the forward scattering amplitude 
for an individual particle ( the cross section da/dfl = \ /(A) | 2 , A= k 2 — \t\ ) and 
JV is a number of particles per unit volume. So, to find the polarization operator IT , 
we have to calculate the forward scattering amplitude normalized in such a way, that it 
reproduces 47r/(0) for an individual particle. Moreover, using the amplitude ( see Sec. 3 
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for its normalization ) T(k\, k 2 ) for an arbitrary momentum transfer A , we obtain for 
x(ki,k 2 ) defined in eq.@ 

X (kx,k 2 ) = T{k 1 ,k 2 )/{2u 2 ). (14) 

From this relation and eq.(^), we can extract c m i(Q,ki) and then, if needed, take into 
account terms with Q ^ omitted at the transition from eq.fllOD to eq.( 



As any other 2x2 matrix, U/(2u) in eq.([13D can be presented in the form of 
n/(2co>) = a + her , where cr are the Pauli matrices. Note that the eigenvalues of 
the matrix a + her are a ± vb 2 • Using a similar representation for the initial density 
matrix p(0) = (I + r]cr)/2 , with 77 being the initial Stokes vector, we find from eq.(|T3|) 

p (L) = ^e~ 7L | - ( cosna + cos/3) + (cosh a — cos/3) - — (g 2 + g 2 ^ + 

+ (gi x g 2 , cr - rj) + (giT?) (g!<r) + (g 2 ry) (g 2 <r) 
- [(g 2 , cr + 77) + (gj x r) , cr)] sin (3 - [(gi , cr + 77) - (g 2 x 77 , cr)] sinh a} , 

(15) 

where 7 = 2Im{a) , a = 2L Im(Vtf) , (3 = 2L Re{Vh^) , g = gx + ig 2 = h/Vh^ . 
The real vectors gi and g 2 satisfy the conditions gig 2 = , g^— g 2 = 1 , since g 2 = 1 . 
As tr p(0) = 1 , the fraction of outgoing photons is given by FRAC = tr p(L) . The 
Stokes vector at the depth L is t](L) = tr cr p(L) / F RAC . In the explicit form: 

i|(L) = | , FRAC = e^ L P , r,\L) = I , (16) 

where 

P = g 2 cosh a — gl cos (3 — (gi , 77) sinh a — (g 2 , 77) sin (3 — 



(cosh a - cos/3)(gi x g 2 , 77) ; 

S = 77 (g 2 cos /? - g 2 cosh aj - [gi - g 2 x 77] sinh a - [g 2 + gi x 77] sin (3 + 
[gi x g 2 + gi (gi , 77) + g 2 (g 2 , 77)] (cosh a - cos/3) . 



(17) 



We emphasize that, according to eq.(^), the polarization degree | tj(L) | is an increasing 
function of depth L except for T7 2 = 1 or P 2 = 1 . 



3 Scattering of photons in a crystal 

We now assume that uo ^> m . Generally speaking, the amplitude of the forward Compton 
scattering off electrons /c(0) may be of comparable size with that due to virtual e + e~ 
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pairs ( VP ) even in GeV energy region ( see discussion in Q ). However, we can neglect 
/c(0) in what follows since it is, first of all, the unit matrix multiplied by some scalar 
factor and, second, real . That is why its contributions to II and IF in eq.(13) are 



cancelled. For the same reason we can omit the real part of the incoherent contribution 
to the forward scattering VP amplitude. Nevertheless, its imaginary part, which does not 
cancel in eq. (|T3|) , will be taken into account. It is the unit matrix multiplied by wW p , 
with W p being the probability per unit length of the e + e~ pair production by a photon. 
Remember that the incoherent probability W p in crystals is smaller than in corresponding 
amorphous media (see ,e.g. ||). So, we must calculate here only the coherent contribution 
to the VP amplitude. 

The properly normalized transition amplitude of the photon with 4-momentum k^ = 
(o»i,ki) and polarization vector into the photon with k 2fl = (c<j 2 ,k 2 ), e 2tl reads 

T(k 1 ,k 2 ) = 2ia J d i x 1 d 4 x 2 Tr[G{x 2 ,x 1 )e 1 e- lklXl G{x u x 2 )e* 2 e tk2X2 ] , (18) 

where a = 1/137 is the fine structure constant. The electron Green function G(xi,x 2 ) 
can be expressed via solutions *ff^ (x) to the Dirac equation in the corresponding external 
field 

iG{x*,x x ) = 9(t 2 -h) £ *W(z 2 ) - eih-h) £ ^-\x 2 )¥-\ Xl ). (19) 

n n 

Recollect that ( see discussion in [|7J ) for uj ^> m the contribution to the amplitude , 
in terms of the non-covariant perturbation theory, is given by the diagram where the 
pair production by the initial photon precedes the annihilation of this pair into the final 
photon . We find, keeping only terms proportional to 9(t 2 — ti) 

T(ki, k 2 ) = 2ia Y / d 4 xid 4 x 2 6(t 2 - t 1 )V nm (x 1 ,e 1 ,k 1 )V+ m (x 2 ,e 2 ,k 2 ) , 

(20) 

V nm (x,e,k) = ¥+\x)ee- ikx ¥-\x) . 



It is really important that in the Fourier-transform of a crystal potential 

U(q) = J dre-^ r U(r) 

only small momenta q |< m are present. For such potentials, the quasiclassical 
operator method can be applied in its standard form to calculate the amplitude T{k\ , k 2 ) . 
Then further transformations in eq.(|20"|) are completely identical to those leading to the 
quasiclassical expression for the probability of e + e~ pair production by a photon in an 
external field ( see Sec. 3 in pi), because V nm in eq.(pO[) is essentially the matrix element 
for this process. Notice that the amplitude T(ki , k 2 ) is the contraction of ei^e^ with 
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the sought tensor T^ v {k\ , k-i) . For the transverse 
components of this tensor we obtain from eq. (pC|) 



with respect to i>\ = \ci/qj\ ) 



rpl] 



i a 



(27T) 



dr exp(— ?Ar) J dr I de 

11 



e(cj — e) 



rf P± (0)e-' A B<- 



.4 



cum 



2e u-e) 



T + 



r/2 



-r/2 



—pi 

m 2 



5 



m 2 + 



v Pl±P2±, 



2s 



to 



CO 



Pl 



.pL. 



Pi±Pl± 



(21) 

where pi^ = p(=Fr/2) being the momentum of an electron on a classical trajectory at 
the corresponding time, A = k2 — ki . When the external field vanishes, the amplitude 
T must vanish as well. It assumes a subtraction which really affects only the term 
proportional to m 2 
below. 



in B lJ . This subtraction will be performed in the explicit form 



Generally speaking, for further calculations we have to know the dependence of p±(t) 
on time ,i.e. the corresponding mechanical problem should be solved. For arbitrary po- 
tential it cannot be done analytically. Fortunately, it is sufficient to know the alteration 
of the quantity (t) during the process formation time Tf . In eq. (pj~D it is the charac- 
teristic size of the variable r contributing to the integral. In the rectilinear trajectory 
approximation (RTA), we substitute into equations of motion r + Vit instead of the 
exact solution r(t) . We refer for details to stating here that in crystals the amplitude 
T can be calculated by the use of RTA for arbitrary photon energy ( oj ^> m ) and crystal 
orientation. 

Representing the crystal potential U(r) as a sum over vectors q of the reciprocal 
lattice 

U(v) = £ G(q) , 



we find within RTA for the transverse momentum 



P±(t) = P±(0) + S(t) = p x (0) - ]TG(q)e»^ - 1 



(22) 



where q\\ = (t'lq) . Substituting eq. (|2~2"[) into eq.(pl|) , we can take the Gaussian- type 
integral over p_i_(0) . Performing the subtraction at vanishing external field and going 
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over to the variable y — 1 — 2e/u , we find 

I* = ^-W!) / drexp(-zAr) f j^J y 

-i y o 



exp H/3r(l+V(r))] 



5^ 



G(q)G(qQ _ (q±q' ± ) e i( q+q ', r ) 



q,q 



y 2 a\ a\ — a l 2 a{ + 5 l [ 



sin{q\\ + g[|)r/2 sinq\\r/2 singer/ 2 



aia 2 ) + m 



(?|| + 9f|W2 gyr/2 

^■0 



9f|T/2 



a li2 = 2 G(q) • ^ 



,«qr 



sinq\\r/2 
W/2 



P 



T- 1p 

get 



2m 2 

u(l-y 2 ) 



(23) 

This expression is still rather complicated for the numerical calculations. It would be 
drastically simplified if we neglect the function ip{ T ) as compared to 1 in the phase of 
eq.(^). We suppose that the angle of incidence $0 with respect to some major crystal 
axis $0^1; since precisely for such angles the strengthening of electromagnetic effects 
happens in crystals as compared with amorphous media. If v>\ is far in azimuth from any 
major crystal plane, then | q\\ | ~| q_L^o I an d ipoff( T ) < (yo/ m /$o) 2 ■ in this estimate 
V denotes a typical magnitude of G(q) being of the order of the axial potential well 
depth. So, in this case we can omit ip(r) in the phase of eq.(|23|) for $ ^> Vo/m. If 
now vi is aligned on some crystal plane, there is a subset of q (q') for which q\\ 
are extremely small or vanish. Those q are perpendicular to the plane. Expanding 
in g|| and qi , we obtain for the corresponding contribution to the phase ip p iij) = 
(r / m) 2 (dU p i(x) / dx) 2 / 12 . Here x is the distance from the plane, U p i(x) being the 



planar potential. As long as iP(t) 
which turns into ippiij) <C k 2 



< 1 , we can substitute r ~ /T 1 into our estimate 
The magnitude of the strong field parameter 

1.32 • 10 ie V /cm 



k(x) 



E„ 



P i 



x 



E 



CO 



m 



m 

±JQ — 

e 



can be estimated using the Table 1 of [§], where, in particularly, averaged over x values of 
(E p i(x)/Eq) 2 are presented for (110) plane of several crystals. As a result, for commonly 
used crystals, ip p i(r) is sufficiently small when u is less than several TeV . Assuming 
that uj and $0 satisfy the formulated conditions, we, finally, rule out ip(r) from the 
phase of eq.(|2"B|) . 



Now the integration over r in eq.([23|) can be performed 

J dr exp [i(q + q - A , r)] = (2vr) 3 <5(q + q - A) 

As a result, the tensor T u acquires the form 

TV(h,k 2 ) = (2Tt) 3 ^2^(Q,k 1 )5(k 2 - h - Q), 
Q 
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where Q = q + q' . Using this presentation of T^ikx^k-z) and eq. (|7|) ,eq. (|I~4"D we find 
that c u (Q , ki) in eq.(^) are c y (Q , ki) = 7r lJ '(Q , fci)/ (4iru) 2 ) . Here we are interested only 
in 7r y (0 , fci) being just the polarization operator II which determines a development of 
the photon density matrix at large thicknesses according to eq . ([13|) . Taking elementary 
integrals over r and y , we obtain from eq.(|23|) 



IT 



auo 



m . E {qi^ KM + <Wl + Qqitf - ) KM - diOx)] 



where /i = 2m 2 /(w | |) and functions are 



(24) 



di(x) = x 2 



1 + x — -x 2 ) g(y/l — ,r) 



[l + x)VT^ 



9(1 -x) 



,l,(.r) ( l g (y/T^) + y/T^^j 0(1 -x) , #(x)=ln^±|_ 



d 4 (a:) = — f f ?»(V1 + x) - v/TT 

7T k V I 



X 



+ 



\ 2 

g(Vl-x) + Vl-x) - 



TlX 

T 



x arctan 



0(1 -x)\ 



Vx^l 0(x-l)+ 



2x 2 r 1 / x — 1 

d 2 (x) = d A (x) + {-g(VT+x) ( — — p(Vl + a;) + Vl + x ) - 3- 



^g(y/l-x) 



7T ^2 
x + 1 



g(VTTx) - Vl-x ) + 1 



7T 2 (l + x) 



0(1 -x)+ 



+ 



arctan 



(x + 1) arctan • 



1 



+ - 1 



vx — T \ - i 

where #(x) = 1 for x > and vanishes for x < . 



0(x - 1)} 



(25) 



The photon energy u enters eq.(^) in the combination /i except of the factor uo 2 
in front of the sum. This sum is almost independent of u near its maximum since the 
latter is given by // ~ 1. As a result, the u -dependence of the quantities a and (3 
in eq.fllSD and eq.([L7|) describing the behavior of the photon polarization is reduced to 
the factor uL for optimal orientations. Correspondingly, the optimal thickness L opt 
is roughly proportional to uj^ 1 . As mentioned above, there is a subset of q at perfect 
planar alignment for which q» vanishes ( \i — > oo ) . The contribution of this subset to IT- 7 
, i.e. the yield of the corresponding planar potential obtained from eq.(^J) reproduces 
well known results ( see, e.g. M and literature cited there ) derived within the Born 
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approximation: 

^f(K 2 (x))Ue\ei + 7e 2 e 2 - 

407T v 

where (...) means averaging over the coordinate x , e± is the unit vector perpendicular 
to the plane and e 2 = v\ x ei . What is lost when ip(r) has been omitted in the phase 
of eq. fl23|) , are higher order corrections in crystal potential . 



Using the explicit form of IP- 7 eq. ( |2~3| ) and adding the incoherent yield,we can find the 
quantities presented in eq . (p~5|) , which describes the properties of a photon beam for any 
initial conditions dependent on the crystal thickness. As explained above, the incoherent 
yield is present in eq. (|T5|) only in the absorption coefficient 7 . We can use any basis to 
calculate IP- 7 . Let this basis be formed by two real unit vectors ei and e 2 satisfying 
eie 2 = Viei = v\e 2 = . We should use the same basis to obtain the initial Stokes 
vector r\ . Supposing that V\ is near some axis direction v 3 , we choose e 2 in the plane 
( reaction plane ) containing V\ and u 3 . For the basis chosen, the circular polarization 
degree is given by the magnitude of the second Stokes parameter £(L) = | T] 2 (L) | . If we 
define the angle of incidence $0 as the angle between V\ and v 3 , and (po which is the 
angle between the reaction plane and the (110) plane, then the explicit form of the basis 
vectors reads 

ei = — e^sin^o + e y cos<p , e 2 = v 3 sm$o — (e^cos^o + e y sm ^0) cost? > 

where e x e y = u 3 e x = u 3 e y = and e x is in the (110) plane. In this basis the 
quantities a and b presented in eq.(|T3]) take the form 

. , . au ^ I G(a) I 2 . t , 

A = q 2 [idi(fj) + d 2 (fj)] , Bi = -(eiq ± )(e 2 q ± ) [id 3 (fx) - d 4 (fj)] , ( 26 ) 
B 2 = , B 3 = X - 



e 2 q±) - (eiq±) [id 3 (fi) - d 4 (n)\ . 



Remember that the eigenvalues of the propagation matrix 11/ (2a;) in eq.(^) are a ± -y/b 2 . 
The corresponding eigenvectors e± satisfying the normalization condition | e± | 2 = 1 
are 

ei + reo e 2 - /'e 2 l>i ^ 



1+ I r I 2 Jl+ I r I 2 63 + Vb 2 



In general, the eigenvectors e± are complex. However, they become real when the quan- 
tity r does so. In particularly, when fi is in the symmetry plane of a crystal like (110) 
plane ((po = ), r vanishes and e± coincide with ei 2 . Just such a case ( y?o = tt/2 , V3 
along < 110 > axis) was the only orientation considered in [fj where a quantity accounting 
for the polarization conversion was calculated by the use of dispersion relations. In our 
notation it corresponds to the term in eq.(E¥) proportional to d^(fi) . We were unable 
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to reproduce the factor in eq.(5) of 0, nevertheless, we emphasize that the expression in 
braces of the cited equation coincides with d^fi)/ fi 2 . 

Consider now, as an example, fully linearly polarized (772 = , rj 2 = 1 ) initial photon 
beam. An efficiency of the polarization conversion process is determined not only by £ (L) 
but also by the fraction of surviving photons FRAC (L) . We use the criterion of such an 
efficiency suggested in || 



FOM(L) = W-£(L)^FRAC(L) , 

L being the crystal thickness. Recollect that £(L) and FRAC(L) depend also on u> 
and V\ . For given orientation of a crystal , we still have a free parameter which is the 
angle of the initial polarization vector with respect to the basis vectors e x and e 2 . We 
do not claim here to the final analysis of the polarization conversion process, so that the 
yield for different orientations will be compared at the same initial condition. Namely, 
we set r\\ = — 1 which corresponds to 4> = n/4 with respect to e 2 ( 4> = 37r/4 with 
respect to B\ ). Evidently, this is the best choice for the alignment on the symmetry plane 
of a crystal when the basis vectors are the eigenvectors of the matrix IP- 7 as well. 

In Fig.l the maximum values of the figure of merit ( FOM ) are shown near < 110 > 
axis of a diamond crystal for uo = 100 GeV as a function of the orientation (angles 
$0 and fo ) . For each direction ( $ , ip ) , the figure of merit was calculated first as 
a function of L , then its maximum value was found. Just these maximum values of 
FOM are plotted, so that different directions correspond usually to different optimal 
thicknesses L opt . Owing to the crystal symmetry, there is no need to perform calculations 
for ip beyond the interval chosen, since they will simply reproduce the results already 
obtained for corresponding (p within the interval. For diamond crystal, the largest effect 
is achieved for $ = l-5mrad off the < 110 > axis on the (110 plane (ipo = 0), where 
( see Fig.l ) FOM ~ 3.7. However, in this case the size of L opt ~ 4.7cm seems to be 
too large for practical use. Generally, a typical size of the optimal thickness is of a few 
cm . For the same axis of a silicon crystal ( see Fig. 2 ) , the peak is at $0 = 2.3 mrad and 
tfo = with FOM ~ 2.6 . We emphasize rather narrow angular with of the peak in 
both directions. The optimal thickness in this case is L opt ~ 8 cm . We have performed 
the same kind of calculations for three major axes of diamond, silicon and germanium 
crystals. Comparing types of crystals, we see that the effect is the largest for diamond 
and the smallest for germanium, where FOM ~ 2.0 can be obtained at L opt ~ 2.8 cm. 
What about orientations, the < 110 > axis is the most preferable, while two others give 
a comparable but smaller yield. 

For (po — the position of a peak is determined by the condition fi = 2m 2 /(u \ 
gii I ) = 1 for the smallest non-zero | ?ii | ■ This is connected with the threshold behavior 
of the functions d\^(fi) in eq.(|24]) at \i — 1 . From this condition we obtain for fcc(d) 
structure near < 110 > axis $™ ax = m 2 l c /(Truj) , where l c is a lattice constant. In 
particularly, for Si we have d™ ax (mrad) ~ 229/u>(GeV) . This fact is illustrated by 
Fig. 3 where FOM (upper curves), FRAC (lower curves), and the degree of circular 
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polarization £ are shown for (fo = near < 110 > axis of a 10cm thick silicon crystal 
as functions of u . Three sets of curves ( from the left to the right ) in Fig. 3 correspond 
to the angles of incidence d^mrad) = 2.29 , 2.08 and 1.91 respectively. All curves in 
Fig. 3 have peaks exactly at the positions prescribed by the condition obtained above. For 
a given orientation ( at fixed $ ), rather narrow shape of these peaks does not allow us 
to handle with the same efficiency a photon beam having the wide energy spread. 

The magnitude of FOM(L) being proportional to £(L) diminishes for partially po- 
larized initial photon beam roughly proportionally to | rj | < 1 as compared to the fully 
polarized case ( | V) | = 1 ). Depending on the orientation, more or less noticeable change 
of the polarization degree | T](L) | occurs for | r) | < 1 . It can be seen in Figs. 4, 5 where 
the absolute values of three Stokes parameters and the polarization degree are presented 
as functions of the silicon crystal thickness L at uo = lOOGeV . The calculations were 



carried out using eq. (|I~6|) , eq. (|TT|) , and eq.(|26|) for r/i = —0.5 , 772 = 773 = . The angle 
of incidence $0 = 2.29 mrad is the same for both figures. At (po — ( Fig.4 ), ^(X) 
is small and | f](L) | ( curve "tot" ) practically does not change due to the smallness of 
the parameter a ( see eq. (|T5|) ) within the whole interval of L presented. The fracture 
of the curve (1) means only that i]i(L) changes its sign at L — 21cm . Since the polar- 
ization degree practically conserves and rjz(L) can be neglected for this orientation, we 
could measure the linear polarization f]i{L) to determine the circular polarization 772 (L) 
appeared. However, the situation drastically changes already for a small alteration of 
(po . It is seen in Fig. 5 calculated at (fo = 0.015 . Now the measurement of the linear 
polarization does not help in the determination of the circular one since the polarization 
degree is no more constant. 

In conclusion, formulas derived describe the propagation of hard polarized photons 
through crystals. They are valid in a wide photon energy range for any orientation and 
any crystal type as long as the approximations used are correct. Our calculations show 
that the linear polarization of multi-GeV photons can be converted with an appropriate 
efficiency into the circular one using properly chosen single crystals. However, if we do 
not use theoretical results for some quantities involved ( e.g. for the polarization degree ), 
the only way to determine the circular polarization appeared in a crystal is the direct 
measurement of it. 
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Figure captions 



1. Fig.l The maximum values of FOM ( the figure of merit , see text ) near < 110 > 
axis of a diamond crystal for uo = 100 GeV dependent on the angles of incidence 
i?o and ipo . 

2. Fig. 2 The same as in Fig.l but for a silicon crystal in the neighborhood of the point 
$o = 2.3mrad , ipo — . 

3. Fig. 3 FOM (upper curves), FRAC (lower curves), and the degree of circular 
polarization £ for tpo — near < 110 > axis of a 10cm thick silicon crystal as 
functions of u {GeV) . Three sets of curves ( from the left to the right ) correspond 
to the angles of incidence do(mrad) = 2.29 , 2.08 and 1.91 respectively. 

4. Fig. 4 Absolute values of the Stokes parameters | rji(L) \ (curve 1), | 772 (L) \ 
(curve 2), | rj 3 (L) \ (curve 3), and the polarization degree | rj(L) \ (curve "tot" ) 
depending on the thickness L of a silicon crystal at u — lOOGel^ , 771 = —0.5 , 
V2 = ^3 = , $0 = 2.29 mrad , if o — 0. 

5. Fig. 5 The same as in Fig. 4 but for (p = 0.015 . 
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